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INTRODUCTION 

Bautin approach to the bifurcation theory of limit cycles has been recently gener- 
alized in the framework of complex analysis ([8] [9]). There are now more cases where 
the Bautin ideal is known. A systematic study of the Poincare center focus problem 
via Abel equations entailed several new examples of Bautin ideals ([2], [3]). This article 
deals with Lienard equations which have been used in many applications (cf. [10]). 
Lienard equations play certainly a key role in Hilbert's 16 th problem, as suggested for 
instance in ([19], [20]), because of the topological simplicity of the return mapping . 
Limit cycles encircle the origin and are necessarily contained in the domain of existence 
of the return mapping. This domain of existence may of course not be equal to the 
domain of convergence of the analytic series which defines the return mapping in a 
neighborhood of the origin. Nevertheless it is interesting to produce an estimate of the 



size of this domain of convergence. Many contributions have been previously done to 
Lienard equations and particularly by N.G. Lloyd and his co-workers ([1], [5], [14], [15] 
,[16], [17]). But the approach via the Lyapunoff series which is used in all these ref- 
erences does not entail information on the domain of convergence of the first return 
mapping itself. We develop here a direct computation of the return mapping (indeed 
of its converse) which yields an estimate of the domain of convergence by a majorizing 
series techniques and a recurrency relation for the coefficients which entail the Bautin 
ideal. Application of the complex analysis methods of ([8], [9]) is then straightforward. 
In Lienard case, the use of the Hironaka polynomial division theorem is replaced by a 
very simple and explicit argument. This displays a bound to the number of complex 
limit cycles in a fixed neighborhood of the origin. The maximal number of complex limit 
cycles on this disk equals 2n. Examples produced with perturbation theory (first-order) 
of Hamiltonians by A. Lins Neto, W. de Melo and C.C. Pugh ([12]) show that this 
bound is optimal (recall that the complex bound is at least twice the real bound). To 
end this introduction we mention other recent references to generalized Lienard equa- 
tions^], [7], [11], [13], [18]) to which eventually the same type of analysis should apply. 
The author thanks M. Caubergh for a careful reading of the manuscript. This article 
was completed during a visit to CRM of Montreal. The author is grateful to L. Belair, 
C.Rousseau and D. Schlomiuck for their invitation and to V. Kaloshin and S. Yakovenko 
for useful comments related to the presentation of the results. 

I-COMPUTATION OF THE BAUTIN IDEAL 

The perturbation of the linear part into a focus adds one real limit cycle to the 
number obtained when the linear part is a center. We consider only the case where the 
linear part is a center to simplify the computations. This paragraph is concerned with 
the polynomial vector fields X of following type 
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(1.1) 
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and the associated flow, solution of the system: 



x = y, (1.2a) 

d 

y = -x + Yy^iAv = -x + p(x)y. (1.26) 
i=i 

The equations of the flow yield a second order differential equation classically named 
the Lienard equation. Note that this differential system is sometimes written: 

x = Y-P(x), (1.2c) 

Y = -x, (1.2d) 

which is easily changed into the preceding ones with y = Y—P(x), if P'(x) = —p(x). 
Write (1.2) in polar coordinates (r, 9): 

x = rcos6>, y = rsin6*. (1-3) 

This displays: 

2rr = 2(xx + yy), rr = r 2 p(rcos6>)sin 2 6> (1-4) 

9 = (xy — yx)/(x 2 + y 2 ) = — 1 + sin6'cos6'p(rcos6'). (1.5) 

This yields: 

dr/d9 = rp(r cos9)sin 2 9 /[— 1 + sin0cos#p(rcos#)]. (1.6) 

Bautin's approach is based on the study of solutions of (1.6) r = r(9) so that 
r(0) = ro, given as an expansion: 



r = r + v 2 (9)r 2 + ... + v k {9)r k + . 
Comparison between (1.6) and (1.7) yields: 
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(1.7) 



^2v' k (9)r k [-l + sin6>cos#p(rcos6>)] 

k>i 

+ ^2[kv k (6)p(rcos6)sm 2 6]r k = 0. (1.8) 
k>i 

This displays the following recurrency relation on the coefficients Vk(9): 
d 

v ' k ($) = Aj(cos0)'sin0[cos0v£_,(0) + (k- l)sm6v k -i(6)}. (1.9) 
i=i 

Let / be an ideal of R[\i, AJ. It is convenient to denote Vk(0) G I to mean that 
for all values of 9, Vk(0) is a polynomial in the parameters (Ai, Ad) which belongs to 
the ideal I. In the following Ck denotes a sequence of non-zero numbers (independent 
of the parameters) which will be defined inductively. Choose the initial conditions for 
the recurrency relation (1.9) as: 

Vl (9) = l. (1.10) 

Then (1.9) yields: 

v ' 2 {9) = Aicos#sin 2 #, thus : v 2 (9) = Ai«4(sin#), u4(0). (1.11a) 

This entails: 

v 2 (2n) = 0. (1.116) 

The first coefficient which really matters for the Bautin ideal is the next one. The 
recurrency relation (1.9) displays: 

v ' 3 (9) = X 1 cos9sm9[cos9v' 2 + 2sm9v 2 ] + A 2 cos 2 #sin 2 # = 0. (1.12) 
This yields: 



v 3 (9) = f 3 (9)+X 1 w 1 3 (sm9), 
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(1.13a) 



with: 

«4(0) = 0, (1.136) 

and 

MB) G (A 2 ), (1.13c) 

and 

/27T 
cos 2 #sin 2 #d# ^ 0. (1.14) 
o 

We prove now by induction the following: 
Lemma 1.1 

Let ko be the maximal integer so that 2ko < k — 1 and k\ be the maximal integer 
so that 2k\ + 2 < k. The coefficient i> k(9) displays the decomposition: 

fci 

i>fc(0) = A(^) + ^A 2j+1 ^ +1 (sin^A), (1.17a) 
i=o 

where: 

f k (9)e(X 2 ,...,X 2ko ), (1.176) 
W 2j+1 (0,A) = 0. (1.17c) 

Proof 

This is certainly true for V2(9) and vs(9). Assume this is so inductively. The 
recurrency relation (1.9) shows that the term wffj^ 1 (sin9) contributes either (for / even) 
to an element of the ideal (A 2 , A2fc ) or (for I = 2h + 1 odd) produces: 

A 2 / l+ i(cos6') 2/l+1 sin6'[cos 2 6'«; 2 !+ 1 (sin6') / + sutfwj^^sinfl)], (1.18) 
which (once integrated against 9 contributes to A 2 / l +iw 2/l+1 (sin6')). 
An immediate consequence of the proposition is the: 
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Lemma 1.2 

The coefficients Vk(2ir) belong to the ideal A2, A 2 fc . 

This last result can be improved for the coefficients Vk(2n) of odd order k. Denote 
n = [d/2] (integer part of d/2). 

Lemma 1.3 

For all odd values of k = 2/c + 1 (ko = 1, n), the coefficient 
v 2k +i(2ir) is such that: 

V2k +i(2n)e(X 2 , A 2fco _ 2 ) + C2fc 1 +iA 2 fco> (1.19) 

with: 

c 2ko +i= / (cos^) 2fco (sin6») 2 ^^0. (1.20) 
Jo 

Proof 

In the recurrency relation (1.9), the only term which contributes to A 2 fc is: 

A 2fco (cos#) 2fcH sin#[sin#], (1.21) 

which yields (1.19). 

The recurrency relation above yields finally: 
Theorem 1.3 

The Bautin ideal of the Lienard vector fields (1.1), defined as the ideal generated 
by the coefficients of the return mapping (or of its inverse) Vk{2n) is equal to the ideal 
(A 2 , A 2n ). The Bautin index B defined as the first integer k so that the increasing 
sequence of ideals generated by the k first coefficients Vk(2it) becomes stationary is equal 
to 2n + 1. 

The recurrency relation entails as well the following result: 
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Theorem 1.4 



For all values of A, there is a neighborhood of the origin on which the number of 
real limit cycles is less than n — 1 (n if we consider the perturbation of a focus). 

Proof 

This is a consequence of the classical Bautin's argument. Collecting the terms of 
the first return mapping, we write the equation for the real limit cycles as: 

^ 3 (27r)r 3 (l + ...) + ^ 5 (27r)r 5 (l + ...) + ... + v 2n+1 (2n)r 2n+1 (1 + ...) = 0. (1.22) 

Successive applications of Rolle's lemma show that the number of real positive zeroes 
of (1.22)is less than n — 1. This shows that the maximal number of limit cycles which 
can bifurcate when the linear part is of focus type is n (because it adds up one term 
in (1.22). Note that this bound is in agreement with the bound foreseen in A. Lins 
neto-W. de Melo-C.C. Pugh conjecture, (cf [12]). Of course this computation does not 
entail any control of the size of the domain on which the number of limit cycles is less 
than n in terms of the coefficients of the perturbation. This is the reason to develop 
further analysis based on the complexification of the return mapping. 

II- Estimates of the radius of convergence of the first return mapping 

Let f\ (x) = J2cik(\)x k be an analytic series in x with polynomial coefficients in 
the parameters A = (Ai, A^). Denote | | (norm of the polynomial a^) as the sum 
of the absolute value of the coefficients and | A | = | Ai | +...+ | Ad |. Recall now the 
following: 

Definition II. 1 

The series fx is called an A -series if the following two conditions are satisfied: 
There are positive constants Ki, K 2 , K$, K 4 such that: 
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deg(a k ) < K x k + K 2 , 



(2.2a) 



I a k \<K 3 Kl (2.26) 

Ao-series form a subring of the ring of formal power series in x with polynomial 
coefficients in A. All the usual analytic operations, like substitution to a given analytic 
function, composition, inversion,... transform Ao-series into themselves. Ao-series have 
been precisely introduced (in the subject) by M. Briskin and Y. Yomdin ([4]). 

In the following, we also denote fx the complex analytic function defined for all 
A G C D on a disc D(0, R) by the A -series. 
In Lienard case, the following holds: 

Proposition II. 2 

The series (1.7) is for all values of 9 a A -series with: 

K x = 1,K 2 = -1, K 3 = tt/2, K 4 = 2. (2.4) 

Proof. 

The inverse series of (1.7), which is indeed the Bautin series: 

r = r + w 2 {6)rl + ... + w k {6)r* + ... (2.5) 

solves the differential equation (1.6). We use Siegel majorizing series techniques. 
A rough majorizing series (both in A and ro) of (1.6) is provided by the solution M of 
the differential equation: 



dM/dO - [M\ A |M/(1 - M)]/[l - | A \M/(1 - M)], (2.6) 
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which displays the expansion: 



M = r + M 2 (9)rl + ... + M k (9)r k + ... (2.7) 

The coefficients M k (6>) are inductively defined by a recurrency relation which is of 
following type: 

M' k {9) = S k [M 2 (0), Mfc_i(#)], k > 2, (2.8) 

where S k is a polynomial with positive coefficients. 
Denote 

M k = Max[M fc (#), # E [0, 2?r]]. (2.9) 
The equation (2.8) yields the following inequality: 

M k < (27r)S k [M 2 , M fc _i], k>2. (2.10) 

This displays a new majorizing series for (2.5) W(r ) solution of the algebraic 
equation: 

W(r ) - r = (2tt)| A \W(r ) 2 /[l - (1 + | A \)W(r )}. (2.11a) 

The algebraic equation (2.11a) has a unique analytic solution W(tq) which is tan- 
gent to ro for small values of tq. This equation yields tq in terms of W{tq): 

ro = W(r ) - (2tt)| A \W(r ) 2 /[l - (1 + | A |)W(r )]. (2.116) 

Elementary considerations on majorizing series show that the converse of a ma- 
jorizing series provides a majorizing series for the converse. 

This immediately provides an estimate of the radius of convergence of (1.7): 



fl(A) = l/[l+| A|], 
and the proof of the proposition II-2. 
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(2.12) 



The estimate of the radius of convergence can be improved (for some values of the 
parameters) with an elementary scaling argument in the case of Lienard's equations. 
The proof given above works actually for an arbitrary polynomial perturbation of the 
rotation flow (with some changes of notations). 

Theorem II.3 

The return mapping of the vector field X converges at least on the disc _D(0, R(p)) 
of radius R(p) = p/2 where p is the unique positive real number such that: 

p d \ A d |+... + p|Ai |=1. (2.13) 

This last estimation improves (2.12) for small | A | but not for large ones. 
Proof: 

Change of coordinates (x,y) into (px,py) transforms 

x = y 

y = - X + p ( x )y^ 

into 

x = y 
d 

y = - x + ^(A iP >V- (2.146) 

i=i 

This means that any result obtained on the disc D(0,R) for the equation (2.14a) 
is valid on the disc _D(0, pR) for the new equation A^ i— > A- = p l \. 

Ill-Estimates of the number of complex limit cycles on a fixed neighbor- 
hood of the origin. 
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(2.14a) 



The generalization of Bautin's approach to complex limit cycles has been done in 
([8], [9]). It yields quite explicit results in the case of homogeneous perturbations of the 
linear part. The same type of results can be displayed in the case of Lienard equations 
and this is done in this paragraph. 

It is quite interesting to note now that in comparison to the "polynomial Hironaka 
division theorem" that we need to use in the general situation, we have to check a very 
easy proposition. 

In Lienard case, the Hironaka basis is obviously A2, ■ ■■\2n an d it yields the following: 
Proposition III.l 

Let / be the Bautin ideal. Let /(A) be an element of I of degree k, there is a 
decomposition: 

n 

/(A) = X>(A)A M , (3.1a) 

i=i 

with 

degifa) < deg(f) -l = k-l, (3.16) 

and 

I <t>i \<\ / I ■ (3-lc) 

Proof. 

Consider /(A), collect in front of A 2 all the monomials which contain A 2 . The 
difference still belongs to the ideal I. Then repeat the process with A4,...The two 
majorations are obvious. 

Following the techniques of ([8], [9]), we prove now that the first return mapping 
belongs to a Bernstein class. It is convenient at this point to change notations and 
write: 
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/(A, r) = r + v 2 (27i)r 2 + ... + v k (27v)r k + (3.2) 

for the (converse of the) first return mapping and /fc(A) = ffc(27r) for its coefficients. 
Indeed, we consider the analytic extension of the series to the complex domain and 
denote now x as the complex variable in place of the real variable r. 

/(A, x) = x + f 2 (X)x 2 + ... + f k (X)x k + ... (3.3) 

Proposition III-2 

The analytic series (3.3) belongs to the Bernstein class B 2 B R ^ with: 

R = K 4 (l+ | A |), (3.4a) 
c = nKiKj n - 1 /c% n+1 . (3.46) 

Proof. 

Write first the condition for f\(x) to be an A -series as follows: 

deg[f k ] <k-l, 
| f k \<K Z K\. 

Write next the decomposition: 

n 

/ fc (A) = J> M (A)A 2i ]. (3.5) 
i=i 

This yields: 



•>k,iW \<\ fk I (1+ I A |) fc " 2 < K 3 Kt\l+ I A \)] k ~ 2 . (3.6) 



This entails: 
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I f k (X < K 3 Kt\l+ I A |)] fc - 2 (Max l=1 ,..., n | X 2l |). (3.7) 

A more careful analysis of the recurrency relation (cf. 1.19, 1.20) shows that equa- 
tion (3.5) entails a nxn matrix relation between the vectors f2j+i,j = l,---,n and 
\2i, i = 1, n of the form: 

n 

f2j + l = + Aji]\2i, (3.8) 

i=l 

where the matrices C and A are respectively diagonal, with non-zero coefficients 
C2i+i defined in (1.20), and nilpotent (upper-triangular). Inverting the matrix relation 
(3.8) yields: 

I A 2 , \<\(C + A)- 1 | Maxi = i,... in | / 2i +i(A) | . (3.9) 
This is completed with the inequalities: 

| (C + A)- 1 \<\ C' 1 || (l + C^A)' 1 |, (3.10a) 

<| C" 1 | [1+ | C~ X A | +... + (| C~ X A \y~ x \. (3.106) 
The coefficients of the diagonal matrix C entail (cf. 1.20): 



C~ L |= l/c 2n+1 , (3.11) 



This yields: 



| (C + A)' 1 |< (l/c 2n+ i) n nK 3 Kl n (l+ | A |)] 2n+1 . (3.12) 
Equations (3.12), (3.9), (3.7) entail now: 



/fc(A) | R k < (l/c2 n+ i) n nKl[K4l+ | A |)] 2n - 1 Max Jfc=2i ... iB (| fk(X) |), (3.13a) 
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I fk(X) I R h < (l/c 2n+1 ) n nKl[K 4 (l+ | A |)] 4 "Max J=2 ,..., B (| fj(X) \)W . (3.136) 

This means that the analytic series (3.3) belongs to the Bernstein class B^ B R ^ 
(cf. [8], [9]) with: 

c = nKl[K 4 (l+\ X\)] 4n /4 n+1 . 

The proposition III. 2 and the results of ([8], [9]) now imply the following: 
Theorem III.3 

The number of zeros of f\(x) in the disc D(0, R') is less than B — 1 = 2n with 

R' = c^ n+1 /[2^nK![K,(l+ | A |)] 4 - +1 , (3.14a) 

R' = C2n+i/[2 10n - 1 n 7 r 2 (l+ | A |) 4n+1 . (3.14a) 

For small values of | A |, this estimate can be improved with the same scaling 
argument as used in the previous paragraph as follows: 

Theorem III.4 

The vector field X has less than 2n complex limit cycles on the disc _D(0, R'{p)) of 
radius 

R'{p)=pc n 2n+1 /[K 2 n2 l ^l (3.15a) 
where p is the unique positive real number such that: 

p d | A d | +... + p | Ai |= 1. (3.156) 
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